Transport Properties Calculation in the Superconducting State for a Quasi-Twodimensional System 
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We performed a self-consistent calculation of the transport properties of a d-wave superconductor. We used 
for calculations the T-matrix approximation. The coresponding equations were evaluated numerically directly 
on the real frequecy axis. We studied the ab-plane charge dynamics in the coherent limit. For the c-axis charge 
dynamics, we considered both, the coherent and the incoherent limit. We also have calculated the penetration 
depth in this model. 
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I. INTRODUCTION 

Although, a large number of experimental and theoreti- 
cal investigations indicate that two dimensionality of the nor- 
mal and superconducting state is one of the key factors in 
high temperature superconductivity [|l|,^, it remains to exam- 
ine whether the cuprates are two-dimensional {2D) metals or 
three-dimensional (30) metals with strong anisotropy. The 
high temperature superconductors show remarkable deviation 
from Fermi liquid behavior in the normal state, such as the 
appearance of the pseudogap phenomena and their related is- 
sues. The pseudogap phenomena mean the suppression of the 
low frequency spectral weight without any long range order. 
There are enormous studies from both experimental and 
theoretical point of view in order to explain these phenom- 
ena. The direct measurements of electronic spectrum such as 
ARPES have indicated the similarity between the pseu- 
dogap and superconducting gap while using Intrinsic Tunnel- 
ing Spectroscopy Q, it was found that the pseudogap, is co- 
existing with the superconducting gap, indicating a different 
nature of the two phenomena. The crucial test for supercon- 
ducting origin of the gaps is their magnetic field dependencies. 
Magnetic field is a strong depairing factor and destroys su- 
perconductivity when the field exceeds the upper critical field 

Hc2- 

The opening of the pseudogap has drastic effect on the 
physical properties of the high Tc cuprates. It is found that 
associated with the pseudogapthe in-plane resistivity deviates 
from the T-linear behavior and the T coefficient of the 
c-axis resistivity changes sign, signaling semiconductor like 
behavior [||]. A remarkable point of the pseudogap is that 
it's structure in momentum space is the same as the supercon- 
ducting dj.2_y2 symmetry with continuous evolution through 
Tc This implies that the pseudogap phenomena have 
close connection to the superconducting fluctuations [p^, and 
strongly suggested that the pseudogap is a precursor of the su- 
perconductivity [|ll[| . The main difference between the super- 
conducting gap and the pseudogap is that the gap function is 
caused by the superconducting order, while the pseudogap is 
caused by the superconducting fluctuations. The optical con- 
ductivity is one of the quantities which most evidently display 
the anisotropy of the system. Measurements of the a6-plane 
conductivity suggest that the conductance in plane is coherent 
both in underdoped and overdoped regimes while the c-axis 



conductivity changes in character from a coherent behavior 
in overdoped regime to an incoherent behavior in the under- 
doped regime JT^]. Theoretical investigations of the a6-plane 
and c-axis conductivity both in the normal and superconduct- 
ing phase can be found in [pj[|. 

In this paper we extend the self-consistent T-matrix calcu- 
lation to the superconducting state. We consider the effects of 
the superconducting fluctuations on the electronic state which 
are the origin of the pseudogap. In the superconducting state 
the T-matrix approximation includes both the amplitude and 
the phase mode of the superconductor order parameter. We 
investigate the transport properties which show the pseudogap 
phenomena. We calculate their behavior in the superconduct- 
ing state. The normal state behavior was calculated in Jl4|]. 
This paper is constructed as follows. In section II we give the 
Hamiltonian and explain the theoretical framework. In section 
III we present the results for the spectral function A(k, uj) and 
the results for the transport properties studied in this paper. In 
section IV we present the conclusions. 



II. THEORETICAL FRAMEWORK 

A. Model Hamiltonian 

In this section we present the theoretical framework of this 
paper. We consider a microscopic model, which incorporates 
both strong electron fluctuation in the CuO planes and a weak 
interlay er coupling. The hopping between layers is included 
in the following Hamiltonian: 
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(1) 



where tc is the hopping amplitude between layers and 
ci,i.ij{cfi 0-) is the annihilation (creation) operator for an elec- 
tron within the planar site i, with spin a, in layer I. Within 
each layer we consider the following two-dimensional model 
Hamiltonian which has dx2-y2 symmetry superconducting 
ground state: 



(2) 



- XI ^k,k'C^k,T'^'.k+q,TC,';k',iC;,k'-q,i 
k,k',c7 
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is the annihilation (creation) operator for 
omentum k, and spii 
electron dispersion relation is given by: 



where c, . (c^ 

an electron with momentum k, and spin a in the layer I. The 



~2t (cos fcr + COS h 



■V) 



At' cos kx cos ky — fj, 



(3) 



where t and t' are the nearest and the next-nearest neighbors 
hopping amplitudes and /i is the chemical potential. For now 
on we will consider t equal to unity and t' — — 0.5i. The pair 
interaction is written as: 



k,k' 



= Vf, 



k,k' 



where 



/k = cos kx — cos ky 



(5) 



is the dx2^y2 wave factor In equation (^ F is negative. The 
important character of high Tc cuprates is the momentum de- 
pendence of the interlayer hopping matrix element ic(k). The 
transfer matrix tc(k) obtained by the band calculation 1 15| is 
expressed as: 



tc(k) — (coskx — coskyY 



(6) 



The large magnitude of the resistivity anisotropy in the normal 
state reflects that the c-axis mean free path is shorter than the 
interlayer distance, and the carriers are tightly confined to the 
CuO planes, and also is the evidence of the incoherent charge 
dynamics in the c-axis direction. 



B. Self-Consistent T-Matrix Approximation 

In this paper we focused on the calculation of the Green 
function directly on the real frequency axis [|l6|| in order to 
avoid the difficulties of controlling the accuracy of the cal- 
culations used in the numerical analytical continuations from 
the imaginary to the real axis by Fade algorithm [[ij]. In the 
previous paper [|l^ we presented the self-consistent equations 
set, which must be solved in order to obtain the Green func- 
tion in the normal state. In this section we will present the 
equations only for the superconducting state. In the supercon- 
ducting state, the self-consistent T-matrix approximation is a 
conserving approximation in the sense of Bay m and Kadanoff 

Si- 

For a continuum model with local interactions, the corre- 
sponding equations have been derived by Haussmann Jl^ ] and 
the analogous equations for the Hubbard model by Pedersen 
et. al. [^o|]. The equations were solved in pi] ] and the ther- 
modynamics of a superconductor was also studied. A similar 
model was also studied in [ ^2[ ] . 

Here we extend the self-consistent T-matrix calculation to 
the superconducting state. We carry out a self-consistent cal- 
culation for the specti-al functions: 



A{k,uj) = ImG{k,Lu) 

TT 



(7) 



B{k,uj) = --ImF{k,uj) 

TT 

where the Green function G(k, oj) is given by: 

LO + ek + S(-k, -uj) 



(8) 



(w - Ek - S(k, w)) (w + ek + S(-k, -ij)) - Al 

(9) 

and the anomalous Green function F{k, uj) is given by: 



(4) F (k, uj) = 



(co - ek - S(k, uj)) (c^ + ek + S(-k, -u;)) ~ 

(10) 



where E(k, w) is the retarded self-energy and Ak is the order 
parameter in the superconducting state. We choose for the 
gap function Ak a (ij.2_j,2 -symmetry form Ak — A/k. The 
imaginary part of the self-energy can be expressed as: 

/mE(k, uj)^ flY^ I doj' [f{uj') + n{u; + uj')] x (1 1) 

k' 

A{k',uj')ImTii{k + k',LU + Lu') 

where f{uj) and n{uj) are the Fermi-Dirac and Bose-Einstein 
distributions. The real part of the retarded self-energy can be 
calculated using the Kramers-Kronig relation: 



/duj 
— 
TT 



duj' ImJ:{k',uj') 



(12) 



where p.v. represent the principal value of the integral. We 
have to mention that only the diagonal part of the T-matrix 
is taken into account when we calculated the self-energy. Ig- 
noring the off-diagonal part of the T-matrix means that the 
anomalous self-energy is not considered. The T-matrix is 
given by the following relation: 



T{(i,uj)^V[l + VU{(i,uj)] ' 
where n(q, w) is a 2 x 2 matrix given by: 

' \L*{q,uj) K{-q,-uj) 



(13) 



(14) 



and: 



The imaginary part of K{q, uj) and L(q, uj) can be calculated 
as follows: 

ImK{(i,uj) - v^E/ d^'fi*9h^ x (15) 
A{k,uj')A{q-k,uj~uj') 

/uj' 
duj'fltgh— X (16) 

K 

B{k,uj')B{(i-k,uj~uj') 

The real part of the K{q, uj) and £(q, uj) are also calculated 
using Kramers-Kronig relation. The chemical potential is de- 
termined by fixing the carriers density n through the relation: 
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^2Y, f dLoA(k,Lu)f{Lu) (17) 



We will consider in the following, the hole doping 6 ^ 1 — n. 
We also have the sum rule which must be satisfied; 



(18) 



The gap A can be calculated self-consistently using the gap 
equation which in our notations can be written as: 

A = VJ2 f dojf^B{k,u;)fiu;) (19) 

Equations ([l7[)-(|l9|) must be satisfied at each step in the self- 
consistent calculation in order to obtain the new chemical po- 
tential and the gap function. We have to mention at this point 
that we did not take into account the anomalous self-energy at 
this level of calculations. The gap function can be calculated 
from the equation: 1 + VK{0, 0) - VL{0, 0) = 0, which is 
equivalent to the gap equation and is realized in the supercon- 
ducting state. We determined the gap Ak using both methods 
and found similar behavior as function of T. 

For this calculation we used a Brillouin zone divided into 
32 X 32 lattice and the frequency integration was done over 
1024 points. All the correlations and the convolutions ap- 
pearing in Eqs. ([rT|)-(|l^ were done using the FFT algo- 
rithm. Performing the self-consistent calculation we obtain 
the spectral functions A(k, w) and -B(k, uj). Using this spec- 
tral functions we can calculate different physical characteris- 
tics of cuprates. We can also calculate the density of state: 



iVM = ^A(k,c.) 



(20) 



The a6-plane conductivity can be calculated using the follow- 
ing relation 1 13|: 



d LU 

duj' 



E 

k 



dkr 



dk,. 



(21) 



a„,(c.) = -de'- V / — [fiu' + L0)- fiio')] X (23) 

k,k' •' 

A{k,Uj')A{k',LU + Lu') 

We calculate also, the London penetration depth. This quan- 
tity is proportional with the superfluid density. The in-plane 
penetration depth is given by: 



[A{k, uj')A{k, uj + oj') + B(k, uj')B{k, uj + uj')] 

The equation ( pT| ) for the ab-plane conductivity neglects the 
vertex corrections and is coherent in character For the cal- 
culation of c-axis conductivity we consider both the coher- 
ent and incoherent limits. The coherent c-axis conductivity is 
given by the relation: 

aciLu) = -de'- V tlik) f — [fiio' + fico')] x (22) 

[A{k, uj')A{k, uj + u;') + B{k, Lu')B{k, uj + u;')] 

In equations (pT|)-(p2[) d represents the interlayer distance. 
The contribution from the incoherent process has been dis- 
cussed by different authors [p3|]. The contribution from the 
tunneling process can be written as: 



XI ~ d 



Y^jd. 



3£k 



5£k 

dk„ 



F{k,uj)\^ f {uj) 
(24) 

and the c-axis penetration depth is: 

J_ = Me'- V / dujtl{k) |f (k, uj)\' f {uj) (25) 

In the next section we present the results obtained for the 
transport properties calculated in this paper and compare the 
results with other theoretical works and with experimental 
data. 

m. RESULTS 

In this section we present the results obtained for the self- 
consistent T-matrix approximation set of equations, intro- 
duced in the previous section. Through out of our calcula- 
tion we chose for the coupling constant V — —hM.To study 
the superconducting state, we first calculated superconduct- 
ing transition temperature T^.. The superconducting transition 
temperature was determined as the highest temperature where 
the set of equations (0)-(|l^ has a non zero solution for the gap 
function. Near the critical temperature the algoritm becomes 
unstable because of the rapid growth of order parameter which 
was found also in the FLEX calculations in [p4|], caused by 
the depairing effect due to low frequency spin fluctuations. In 
our case we consider the effect of the superconducting fluc- 
tuations. The effect is included in the retarded self-energy, 
which is determined self-consistently. Using T-matrix approx- 
imation we include both the amplitude and phase mode of the 
superconducting fluctuations in the self-energy [ ^5| ] . The self- 
energy corrections are reduced in the superconducting state 
because the fluctuations are reduced. The T-matrix approxi- 
mation gives a unified description for the normal state in the 
pseudogap region of the underdoped cuprates and the super- 
conducting state [p^. The approximation is not accurate near 
the critical temperature Tc because of the strong suppression 
of the depairing effect due to the pseudogap, followed by rapid 
growth of superconducting gap below T^. The order parame- 
ter grows more rapidly than in the BCS model. We chose 
5 — 0.15 as the hole concentration. We found the BCS crit- 
ical temperature Tbcs — 0.37i and in the case of the self- 
consistent calculations Tc = 0.22i. In Fig.l we present the 
order parameter A as function of temperature. 
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Fig. 1. Temperature dependence of the order parameter A. The 
solid hne represents the BCS gap function while the symbols line 
represent the self-consistent gap function. 

After performing the self-consistent calculation we can cal- 
culate the transport properties of the system. The penetration 
depth can be calculated alon g th e afe-plane and c-axis direc- 
tions, using Eqs. ( ^4| ) and (|25[). In Fig. 2 and Fig.3, we 
present the temperature dependence of the penetration depth 
as function of temperature 



We found a similar behavior as in the case of BCS model. 
Due to the rapid increase of the order parameter blow Tc 
we found a rapid decrease of the penetration depth with de- 
creasing temperature in both a&-plane and c-axis directions. 
Recently, the penetration depth was measured as function of 
temperature by different authors for YBCO and Bi — 2212 
using different techniques [|6|]. The magnetic field penetra- 
tion depth A of superconductors is related to the supercon- 
ducting carrier density Ug divided by the effective mass m* 
as 1/A^ oc ns/m* . In a recent paper Uemura |^] explained 
the experimental data found in [^] assuming a microscopic 
phase separation between superfluid and non-superconducting 
fermionic carriers, similar to the superfluid He films. 




T/t 
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Fig. 2. Temperature dependence of the in-plane penetration depth 
in the BCS model (solid line) and the self-consistent calculation 
(symbols line) 



Fig. 3. Temperature dependence of the c-axis penetration depth 
in the BCS model (solid line) and the self-consistent calculation 
(symbols line) 

Recently, it was reported that, the photoemission spectra 
near the Brillouin zone boundary, exhibit unexpected sensitiv- 
ity to the superfluid density p^], supporting the universality 
as function of doping and temperature of the properties of the 
cuprates. We did not consider the vertex corrections for the 
calculation of i^(k, oj). Desired improvements of the present 
approach include the consideration of the vertex corrections 
neglected in the present approach. 

We also analyze the afc-plane and c-axis optical conductiv- 
ity in the superconducting state using the same approximation. 
For the calculation of Oab we consider a coherent nature of the 
interlayer coupling. For the c-axis conductivity we consider 
the coherent and the incoherent coupling. Coherent coupling 
originates from an overlap of electronic wave functions be- 
tween planes, and in-plane momentum is conserved in inter- 
layer hopping. By contrast, for impurity mediated incoherent 
coupling, the in-plane momentum is not conserved. 

The influence of the nature of the interlayer coupling on c- 
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axis conductivity was studied by different authors yl^,?]. The 
a6-plane conductivity is mainly due to the quasiparticles near 
the 'cold points', although the c-axis conductivity is due to the 
quasiparticles near the 'hot points' . Since the superconducting 
^2-2 _ j,2 gap and the pseudogap in the underdoped region above 
Tc are large at 'hot points', the c-axis conductivity reflects the 
pseudogap more clearly then the a6-plane conductivity. 

In Fig A we present the calculated (Tab, based on Eq. (|2l|). 
The behavior in the superconducting state is similar with the 
behavior in the normal state. 
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Fig. 4. Frequency dependence of the afc-plane conductivity for 
different temperatures. The solid line corresponds to T — 0.19t, 
the dashed line to T = O.lt and the dotted line to T = 0.07t. The 
critical temperature is Tc ~ 0.22t. 

The in-plane spectrum is dominated by a Drude peak at 
uj = 0. The intensity of the peak, increases with increasing 
temperature. Below Tc there is no signature of superconduct- 
ing gap seen in Y Ba2Cu^0j-s [p^. Careful study of opti- 
cal conductivity by different methods shows that the width of 
the Drude peak diminished by several orders of magnitude by 
decreasing the temperature, just below Tc [|T]|. In cuprates, 
in contrast with BCS theory, the excitations are electronic, 
and as the gap develops in this excitations, the decrease in the 
scattering take place. The optical properties of the cuprates, 
are those of a clean limit of a BCS superconductor For the 
calculations of the c-axis conductivity we consider both the 
coherent and the incoherent nature of the interlayer coupling. 

In Fig. 5 we present the results for the c-axis optical con- 
ductivity in the case of the coherent nature of the interlayer 
coupling. The approximation assumes the independent elec- 
tron propagation in each layer and is justified for tc <C t. In 
this case, the (Tc spectrum is dominated by a Drude peak at 
uj = 0. This features are characteristics of the rf-wave sym- 
metry, and due to the contribution from the gap node. 




Fig. 5. Frequency dependence of the coherent c-axis con- 
ductivity for different temperatures. The solid line corresponds to 
T = 0.19t, the dashed line to T = O.lt and the dotted line to 
T = Omt. The critical temperature is Tc = 0.22t. 




Fig. 6. Frequency dependence of the incoherent c-axis con- 
ductivity for different temperatures. The solid line corresponds to 
T = 0.19t, the dashed line to T = O.lt and the dotted line to 
T = omt. The critical temperature is Tc = 0.22t. 

In Fig.6, the results for the c-axis conductivity are pre- 
sented in the limit of the incoherent nature of the interlayer 
coupling. In this case the electronic contribution to CTc {^) 
does not form a Drude peak at — |32|. In the super- 
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conducting state, the c-axis optical conductivity is suppressed 
furthermore with decreasing temperature and shows the gap 
structure at low temperatures. We consider both limit, because 
measurements of c-axis conductivity suggest that conductance 
in c-direction is coherent in the overdoped regime [|33|], and 
incoherent in the underdoped regime [ |l2| ] . The coherent 
was calculated using Eq. (g^) ^nd the incoherent Oc using Eq. 
©. 

The results presented in Fig.6 are in agreement with 
the experimental results obtained for the underdoped 
YBa2Cu30T-s- Here we calculate the conductivity by ne- 
glected the vertex corrections. The vertex correction is not 
important, except for the Umklapp scattering in the case of 
electron correlation [§41. 



IV. CONCLUSIONS 

Solving the self-consistent T-matrix equations for the 
model presented in the section 2 we have computed the trans- 
port properties of the system in the superconducting state. 
Calculations were performed for the superconducting order 
parameter of dx2_y2 symmetry. The effect of the supercon- 
ducting fluctuations is included in the self-energy corrections. 
It is clear that the T-Matrix approximation breaks down in the 
limit T Tc since the propagator is dramatically renormal- 
ized at the Fermi level. Our scenario is based on resonance 
scattering [pi]]. The gap function was found to develop more 
rapidly than in the BCS model. We calculated the London 
penetration depth and the c-axis and a5-plane conductivity 
and compared the results with the experimental data. The 
study of experimental dependence of the penetration depth is 
an important problem [ p7| , p8| ] . 

Desired improvements of the present approach include the 
considerations of the vertex corrections neglected in this ap- 
proximation. The vertex corrections have a significant effect 
in the calculation of Hall resistivity [ ^5[ | and magnetoresis- 
tance [|6|] in high-Tc superconductors. Although, the effect 
was suggested to be insignificant to superconductivity [^7|], a 
self-consistent calculation included vertex corrections is de- 
sired. 
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